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The generalized parton distributions of the pion are studied within different light-front approaches for the 
quark-hadron and quark-photon vertices, exploring different kinematical regions in both the valence and non- 
valence sector. Moments of the generalized parton distributions which enter the definition of generalized form 
factors are also compared with recent lattice calculations. 
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1. Introduction 

Generalized parton distributions (GPDs) rep- 
resent a key concept for understanding the hadron 
structure |1|2|3|4|5] , They unify the information 
encoded in electromagnetic (e.m.) form factors 
(FFs) and ordinary parton distributions, sup- 
plementing them with the possibility to access 
new aspects of the hadron structure. In particu- 
lar, the pion GPDs represent an important test 
ground for model calculations aiming to a de- 
tailed description of hadron structure, and this 
explains the wealth of papers devoted to such 
a task [6171819110111112113114115116] . Here we re- 
view the results of Ref. [16] for the calculation of 
the pion GPDs in three relativistic models which 
explore different kinematical regions in a comple- 
mentary way. In particular, the first model is a 
covariant analytic model, based on 4D Ansatze 
for the pion Bethe-Salpeter amplitude (BSA), 
which allows to explore the whole kinematical do- 
main in the valence and non-valence sector. The 
other two models are constrained to either the va- 
lence or non- valence regions. In the non- valence 
region we adopt a model based on a microscop- 



ical vector-meson model dressing for the quark- 
photon vertex and a phenomenological Ansatz for 
the 3D light-front (LF) projection of the pion and 
vector-meson BSAs. Finally, in the valence region 
we discuss a third model constructed within the 
LF relativistic Hamiltonian dynamics. 
After a short introduction about the general for- 
malism for the definition of the pion GPDs, in 
Sect. [2] we present the main features of the three 
models, referring to [16] for a more detailed dis- 
cussion. In the final section, we show the model 
results for the pion GPDs, together with the com- 
parison between our model predictions and recent 
lattice results for the first moments of GPDs en- 
tering the definition of generalized form factors 
(GFFs). 

2. GPDs in covariant and light-Front rela- 
tivistic models 

GPDs are defined as the non- forward (p ^ p') 
matrix elements of light-cone bilocal operators 
separated by a light-like distance, i.e. 



G r = {p',7r ± \O r b,^}, 



(1) 



1 



2 



with 



dz 
An 



o±). 

(2) 



In Eq. j2]), P = (p + p')/2 is the average pion 
momentum and the operator T is a matrix in 
the Dirac space which selects different spin po- 
larizations of the quark fields. In particular, for 
r = 7 + one has the unpolarized quark GPD H q , 
while r = ia +l projects on the transverse polar- 
ization of quarks and defines the chiral-odd GPD 
Erf . Because of Lorentz invariance the GPDs can 
only depend on three kinematical variables, i.e. 
the (average) quark longitudinal momentum frac- 
tion x = k + /P + , the invariant momentum square 
t = A 2 = (p' — p) 2 , and the skewness parameter 
£ = — A + / (2P + ). In addition, there is an implicit 
scale dependence in the definition of GPDs corre- 
sponding to the factorization scale u 2 . The vari- 
able x allows one to single out two kinematical 
regions. The first region corresponds to the va- 
lence contribution and is given by the union of the 
interval x G [— 1, — 1£|] (for an active antiquark) 
and x £ 1] (for an active quark). In the Fock- 
space expansion of the pion state, this region is 
described by matrix elements with the same num- 
ber of partons in the initial and final states. The 
second region corresponds to x € [— 1£|, |£|], and 
is associated with the non-valence contribution 
involving non-diagonal matrix elements between 
parton configurations with An = 2. 

In the forward case, p — p' , both A and £ are 
zero, and H reduces to the usual parton distribu- 
tion function, while Et vanishes for time-reversal 
invariance. 

Moments in the momentum fraction x play an 
important role in the theory of GPDs. Weighting 
Eq. (fTJ) with integer powers of x and integrating 
over x, the quark operator O r reduces to a local 
operator and the corresponding matrix elements 
can be parametrized in terms of GFFs, i.e. 

r+i n 
/ dxx n - 1 H"(x,U) = J2( 2 2iA U(t), 

/ dxx n - 1 EUx,U)^Y,( 2 2l B q T n.2 i (t)- (3) 
•1-1 i-n 



In Eqs. ([3|), the lowest moment n = 1 of the un- 
polarized GPD yields the pion e.m. FF, while 
the Fourier transform of B^, 1 (t) in the impact- 
parameter space determines the dipole-like distor- 
tion of the quark density in the transverse plane 
due to the transverse spin-structure of the quark 
in the pion. The second Mellin moments of unpo- 
larized GPDs can be related to the GFFs of the 
energy-momentum tensor of QCD. 

The starting point of our approach is the Man- 
delstam formula [17j for the quark correlator in 
Eq. I|2[). giving for the u-quark unpolarized GPD 

H u {x,t,t) = -iN c nJ ^L I S(P+x-k+) 

xA(fc- P,p')A(k- P,p) 
xTr ^S(k - P) 7 5 5(fc + jh+S(k - ^) 7 5 

(4) 

where N c — 3 is the number of colors, 1Z = 
2m 2 / f 2 , with J' 77 the pion decay constant, and 
m and S(p) are the mass and the Dirac prop- 
agator of the constituent quark (CQ), respec- 
tively. In Eq. 75A(fc,p) is the pion vertex 
function deduced from a simple effective quark- 
pion Lagrangian [18j . In the following, we will 
explore different approximations to model the 
momentum-dependent part of the vertex func- 
tion. 

In a first covariant analytic model, A is assumed 
to be a symmetric function of the two quark mo- 
menta with the following form 



A(k-P,p) = C 



(k- A/2) 2 -m 2 R + ie 
1 



(P - k) - m\ + te 



(5) 



A different choice, based on the sum instead of 
the product of the two terms in Eq. (O, was 
adopted in Ref. [19] for the calculation of the 
e.m. FF and further discussed in the case of the 
GPDs in Ref. [T^] • However, the product form ([5]) 
provides a more realistic transverse-momentum 
falloff, leading to a more favourable comparison 
with the experimental data for the e.m. FF at 
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high-momentum transfer and also satisfying the 
support conditions for the parton distribution. 
Once the CQ mass is fixed and the constant C 
in Eq. is constrained through the charge nor- 
malization, the only free parameter of the model 
is the regulator mass mn, which is fitted to the 
experimental value of /„•. The projection into 
the valence and non-valence contributions to the 
GPDs is obtained after integration of Eq. (0| over 
the LF energy k~, fully taking into account the 
pole structure of both the Dirac propagators and 
the vertex functions. 

A second covariant model is introduced by fol- 
lowing the approach of Ref. [20] for the calcula- 
tion of the e.m. FFs both in the spacelike and 
timelike region. Starting from the same formal 
expression of Eq. ((4j) for the GPD, we introduce 
the following new ingredients: i) instead of the 
bare quark-photon vertex, 7^, a dressed quark- 
photon vertex r M (A;, A), modeled through a mi- 
croscopical vector meson (VM) dominance ap- 
proach, and ii) phenomenological Ansatze for the 
BSAs in the valence and non- valence regions. An- 
other basic difference with respect to the previ- 
ous analytic model is that only the simple an- 
alytic structure of the Dirac propagators is re- 
tained, i.e. the analytic structure is disregarded 
in the BSAs of both i) the initial and final pion 
and ii) the VM dressing of the quark-photon ver- 
tex. This approximation turns out to be a very 
effective one in the calculation of the e.m. FF 
only in the A^ = frame [21]. which will be 
also adopted for the present calculation of the 
GPDs. In the valence sector, after integrating 
over the LF energy, the resulting momentum- 
dependent part of the 3D BSA of the pion and 
vector mesons are approximated with light-front 
wave functions (LFWFs) which are eigenstates 
of the CQ square mass operator of Ref. [22 . In 
the non- valence region, there is also a pion non- 
valence component describing the emission (ab- 
sorption) of a pion by a quark. Assuming a van- 
ishing pion mass, such a process can be described 
using a constant interaction, with a coupling con- 
stant fixed by the normalization of the pion FF. 
Furthermore, in the limit m K = 0, only the 
pair-production mechanism is contributing to the 
GPD. In this term, for = one has only in- 



stantaneous contributions produced by the stan- 
dard LF decomposition of the propagator ( i.e. 
S(k) = (Jf on + m)/[fe+(fc--A- +ie)]+7+/2fe+). 
In order to model the instantaneous vertex func- 
tions, we put A lst w CA^" H , where the constant 
C is thought to roughly describe the effects of the 
short-range interactions. Indeed, we use the rel- 
ative weight Cvm/Ck as free parameter. 

Finally, a third model calculation is based on 
a light-front Hamiltonian (LFH) approach based 
on a Poincare covariant description of the pion. 
In particular, the rotational covariance is fulfilled 
through the introduction of the Melosh rotations 
and the proper definition of the total intrinsic an- 
gular momentum. Such a model allows us to ex- 
plore only the valence region, and therefore will 
be discussed just for the £ = kinematics. In such 
a frame, the GPDs can be expressed as overlap of 
LFWFs, given by the product of the momentum- 
dependent part of the wave function in the ini- 
tial and final state with a spin-dependent part as 
dictated from the proper Melosh transformations. 
For the momentum-dependent part we adopt a 
gaussian form [23], with the quark mass and the 
gaussian width fitted to the pion charge radius 
and decay constant. 

3. Results 

In Fig. Q] we show the unpolarized GPD H in 
the (x,t) plane at fixed |£| = 1, for the isoscalar 
u + d and isovector u — d quark combination, 
comparing the results for the covariant analytic 
model with a product-form for the BSA with the 
phenomenological BS model with dressed quark- 
photon vertex. The value |£| = 1 corresponds 
to the contribution of the pair-production mech- 
anism in the whole range of x. The general shape 
for the GPDs in the two models show similar fea- 
tures, with the collinear peak at x ~ 1 increasing 
at higher values of As discussed in Ref. [16] ; 
the covariant analytic model exhibits an overall 
agreement also with the LFH model at £ = 0, and 
can be used at any value of x, £, t for interpolating 
between the other two phenomenological models. 
In particular, at the crossing point of the valence 
and non- valence region, x — |£|, the covariant an- 
alytic model predicts a smooth transition, due to 
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H«>(x,lSl=l,t)/F moll (t) H Jt i=i(x,l^l=l,t)/F mon (t) 




Figure 1. Upper left (right) panel: isoscalar (isovector) unpolarized GPD from the covariant analytic 
model with the product-form for the BSA (Eq. (0) at |£| = 1 and = 0. Lower panels: the same as 
in the upper panels for the microscopic model with dressed photon-quark vertex. On the z-axis the ratio 
with respect to F mon = 1/(1 + \t\/m 2 p ) is presented. The figure is adapted from Ref. [IB] . 



-t (GeV/c) 



2 3 
-t (GeV/c) 2 



Figure 2. Left panel: the ratio A^ (t) / A^ (0) as a function of t. Solid line: product-form for the pion 
BSA, Eq. ([5]), and = 0. Dotted line: the same as the solid line, but with = 140 MeV. Dashed 
line: LFH model, with a gaussian pion wave function and the proper Melosh rotations. Shaded area: 
results from lattice QCD [24] . Right panel: the same as the left panel, but for Alp^ 1 (t) j Alp^ (0) . The 
figure is adapted from Ref. [16]. 



the continuity of the model. 

In Fig. [2] we show results for the ratios of 
the GFFs A{f{t)/Aif{Q) and Ai\f {t) / A{f {0) 
which are evolution-scale independent [13j . The 
predictions from the covariant analytic model 
for two different values of the pion mass are 
shown for both GFFs, while the LFH model can 
be only applied for A^WMfo^ ) at £ = 0. 
The dashed band shows recent lattice results de- 
scribed through a monopolc form 1/(1 — £/Mf j), 
as obtained in Ref. |24j . In particular, we used 
M 2j0 = 1.329 ±0.058 GeV and M 2 . 2 = 0.89 ±0.25 
GeV, corresponding to an analysis of the lattice 
data that satisfies the low energy theorem, i.e. 
^2^°(0) = -44^° (0). Our model predictions 
are overall consistent with the lattice results, ex- 
cept for small values of \t\. A better description 
of the low |i| region could be obtained by incor- 
porating in our phcnomenological models inter- 
action terms responsible for the confinement. On 
the other side, the large uncertainties in the lat- 
tice results for A I 2 = ^ > do not allow us to elaborate 
too much on the comparison between the different 
predictions. 

Finally, in Fig. [3] we show the density for trans- 
versely polarized quark in the impact-parameter 
space b±. Such a density is defined as |25|26j 



2m„ 



db\ 



D 



,o(S±). ( 6 ) 



where the GFFs in the impact-parameter space 
are obtained by Fourier transform of Eqs. ([3]). 
In Eq. ([6]), the monopole distribution associ- 
ated to H is distorted by a dipole term pro- 
portional to Et in the case of transversely po- 
larized quark. The results in Fig. [3] correspond 
to the LFH model and exhibit a clear correla- 
tion between quark spin and transverse space. 
The average sideway shift amounts to (b^)" 11 = 
B^ 10 /(2m w AY t0 ) = 0.197 fm. Remarkably, this 
value is of the same strength as the dipole-like 
distortion in the density of transversely polarized 
quarks in an unpolarized nucleon, i.e. = 
Q /(2m N A" ) = 0.209 fm, as obtained in a 
LFH model for the nucleon [57] • These results 
are also supported from recent lattice calcula- 
tions [SUB], giving (6^} u = 0.151(24) fm and 
(b v ± ) u — 0.154(6) fm for the pion and nucleon, 
respectively. 
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Figure 3. Density in the impact-parameter space 
of transversely polarized u quark in as pre- 
dicted from the LFH model. 
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